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. . , A Gauss-Bonnet dark energy model is considered, which is inspired in string/M-theory and takes 

. also into account quantum contributions. Those are introduced from a conformal quantum anomaly. 

' The corresponding solutions for the Hubble rate, H, are studied starting from the Friedmann- 

' Robertson- Walker equation. It is seen that, as a pure effect of the quantum contributions, a new 

solution for H exists in some region, which does not appear in the classical case. The behavior 
of all encountered solutions is studied with care, in particular, the role played by the quantum 
r-{ ' correction term — which depends on the number of matter fields — on the stability of the solutions 

, around its asymptotic value. It is argued that, contrary to what happens in the classical case, 

^ ■ quantum effects remarkably lead to the realization of a de Sitter stage which corresponds to the 

^sO ' inflation/dark energy stages, even for positive values of the fo constant (coupling of the field with 

the Gauss-Bonnet invariant). 

cr. 

, I. INTRODUCTION 

^""^ ' One of the most intriguing questions in today's Physics concerns the nature of the dark energy known to be 
^ present in the Universe we live in. The existence of such energy — with an almost uniform density distribution and a 
substantial negative pressure — which completely dominates all other forms of matter at present time, is inferred from 
, several astronomical observations [T|. In particular, according to recent astrophysical data analysis, this dark energy 
• seems to behave very similarly to a cosmological constant, which is responsible for the accelerating expansion of the 
observable universe. Furthermore, there are strong reasons to believe that answering this question will have much to 
do with the possibility to explain the physics of the very early Universe. 
pC. ' Models of dark energy are abundant. One of the proposed candidates for it is the so-called phantom field, thus 
named because it corresponds to a negative-energy field. The peculiar properties of a phantom scalar (with negative 
kinetic energy) in a space with non-zero cosmological constant have been discussed in an interesting paper by Gibbons 
0. As indicated there, phantom properties bear some similarity with quantum effects An important property of 
the investigation in ^ is that it is easily generalizable to other constant-curvature spaces, such as Anti-de Sitter (AdS) 
?H space. Presently, there is considerable interest in such spaces, coming in particular from the AdS/CFT correspondence. 
According to it, the AdS space may have cosmological relevance [4j, e.g. by increasing the number of particles created 
on a given subspace [5|]. It could also be used to study a cosmological AdS/CFT correspondence fd\: the study of a 
phantom field in AdS space may give us hints about the origin of this field via the dual description. In the supergravity 
formulation one may think of the phantom as of a special renormalization group (RG) flow for scalars in gauged AdS 
supergravity. (Actually, such a RG flow may correspond to an imaginary scalar.) 

Another candidate for dark energy is the tachyon 0, This is an unstable fleld, actually. The interest of models 
exhibitin g a tachyon is motivated by its role in the Dirac-Born-Infeld (DBI) action as a description of the D-brane 
action [ol llOl. [ll| . In spite of the fact that the tachyon represents an unstable field, its role in cosmology is generally 
still considered to be useful, as a source of dark matter [a [11] and, depending on the form of the associated potential 
[T3I [T^ [Tsl [TgI it can actually lead to a period of inflation. On the other hand, it is important to realize 
that a tachyon with negative kinetic energy (yet another type of phantom) can also be introduced i^Sj. In that 
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phantom/tachyon model the thermodynamical parameter w is naturally negative. In this case, the late-time de Sitter 
attractor solution is admissible, and this is one of the main reasons why it can be considered as an interesting model 
for dark energy Moreover, in order to understand the role of the tachyon in cosmology it is necessary to study 
its effects on other backgrounds, as in the case of an Anti-de Sitter background Q. 

Another remarkable proposal is that the origin of the dark energy could in fact be related with the problem of 
the cosmological constant. One of the most interesting approaches to this paradigm is modified gravity. Actually, 
it is not absolutely clear why standard General Relativity should be trusted at large cosmological scales (thus, 
through an enormous range of orders of magnitude). One may assume that, considering minimal modifications, the 
gravitational action contains some additional terms growing slowly in the presence of decreasing curvature (see, e.g., 
m, [13, HH [13, m, m, [Hi and, for a review, see [11]), which could be responsible for the current acceleration. Within 
such scenario one of the most accepted approaches is the model of modified Gauss-Bonnet (GB) gravity. In this 
model additional terms in the gravitational action are introduced by means of a function which depends on the scalar 
curvature, R, and on the Gauss-Bonnet invariant, G [13] ■ It is possible to demonstrate that such models lead to a 
plausible effective cosmological constant, quintessence, or phantom eras. From these results one can conclude that, 
concerning its role as a gravitational alternative for DE, modified GB gravity may become a very strong candidate 

In the present paper we will consider a Gauss-Bonnet model for gravity, together with the first contributions coming 
from quantum effects. A stringy correction will be added to the ordinary General Relativity action, that is, a term 
proportional to the GB invariant G. We will duely take into account the quantum effects coming from matter, and 
their influence on the stability of the de Sitter universe will be carefully discussed. 



II. THE MODEL 



We consider the model consisting of a scalar field, cf), coupled with gravity in a rather non-trivial way and, as 
advanced, we will also take into account quantum effects. 

As a stringy correction, the term proportional to the Gauss-Bonnet invariant, G, is written as: 



G = i?2 - 4i?^,i?'^'^ + R^^.p^R^^'P'' . (1) 



The starting action is given by 



S = j d\^g [^R- ^9^09^0 - V{4>) + /(0)g| . (2) 

Here 7 = ±1. Action ([2]) has been proposed as a stringy dark energy model in [13]. It is interesting to note that one 
can also add to this system R^ and i?^ stringy corrections [H, [1^ (for a general introduction to modified gravity, 
see [13, [131 )■ Moreover, action ^ is able to solve the initial singularity problem [U [13] For the canonical scalar, 
7 = 1 and, at least when the GB term is not included, the scalar behaves as a phantom only when 7 = — 1 [13], 
showing in this case properties similar to those of a quantum field [3]. In analogy with the model in [34| — where also 
a non-trivial coupling of the scalar Lagrangian with some power of the curvature was considered — one may expect 
that a GB coupling term of this kind can be relevant for the explanation of the dark energy dominance nowadays. 
Doing the same as in Ref. .23] , by varying over 0, one obtains 

Q = ^V^-V'{ci,)+f'{cj,)G. (3) 

On the other hand, performing the variation over the metric g^i, (as in [35]), gives 

+ ^g^^R^ + 7 ~ \g^''dpcjydP4^ + \g^^ {-V{cj>) + /(0)G) 

-2/(0)i?i?^'' + 2V''V'^ (/(0)i?) - 25^'^V2 (/(0)i?) 
+8fiq^)R^pR''P - 4VpV^ if{q^)R''n - 4VpV'^ {.f{4>)R''') 

+4V2 U{mn + 45^''VpV. {nmn - 2f{cl))R''f"'^R%,, + 4VpV. {f{cj))R'^P'''') . (4) 
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And using the identities obtained from the Bianchi identity, 

WpW^R^P"" = W'^R"" - ^WWR + R^P^^Rpa - R^pR^P , 
WpWRP" + WpWRP" = i {W^WR + VVPR) - 2RfP'"'Rp^ + iR^^pR^P , 

VpV^i?"'^ = ini?, (5) 

one can rewrite (g]) as 

= -i {-R^'' + Ig^'R^ + 7 (^i^d^<pd^<P - -^g^^dpcj^dP^ + \g^^ {^V{cj>) + /(0)G) 
-2f{<f>)RRP'' + Af{,p)R^^pR-P - 2f{cl>)R'^P''-R%^^ + Af{^)RPP-- Rp, 
+2 (V^ VV(0)) - 25'''' (V2/(0)) - 4 (VpV''/(0)) i?--" - 4 (Vp W(0)) i?^" 
+4 (v2/(</,)) i?'^'' + Ag^^^ [VpVafm R"" - 4 (VpV,/((/))) i?^""-" . (6) 

The above expression is vahd in arbitrary spacetime dimensions. In four dimensions, the terms proportional to /((/>) 
without derivatives are canceled among themselves, and vanish, since the GB invariant is a total derivative in four 
dimensions. 

Quantum effects can be included by taking into account the contribution of the conformal anomaly, T/i, defined as 
follows: 

TA=^h[F+ ^Di?^ + h'G + h"UR , (7) 

where F is the square of a 4d Weyl tensor and G is the Gauss-Bonnet curvature invariant. They are given by 

F - {l/i)R^ -2R,jR'^ + R,,kiR''''\ 

G = R^ - AR,,R'^ + R,,kiR''''' . (8) 

In general, with N scalar, spinor, iVi vector fields, A''2 {— or 1) gravitons, and A'^hd higher-derivative conformal 
scalars, the coefficients h and 6' take the form 

TV + QNi/2 H- 12iVi H- eiliVa - SiVno 
^ 120(47r)2 ' 

_ N + lliVi/2 + Q2Ni + MlliVz - 28jVHD 

" 360(471)2 ■ (9) 

We have that & > and 6' < for usual matter, except for higher derivative conformal scalars. Notice that b" can be 
shifted by a finite renormalization of the local counterterm R^ , thus b" can be arbitrary. 

In terms of the corresponding energy density [36, 37], pA, and of the pressure density, pA, Ta is given by Ta = 
—pA + ipA- Using the energy conservation law in the Friedmann- Robertson- Walker (FRW) universe, 

PA + iH{pA+PA)^0, (10) 

we can eliminate pA, as 

Ta^~Apa- pa/H . (11) 

This yields the following expression for pA'- 

^ j Ata^HTA 

^ j dta*H -12bH^ + 24b'{-H^ + H^H + H^) 

-{Ab + Qb"){H + 1HH + ^11^ + 12H'^H^ . (12) 



PA = 

= I Ata^H 
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Now, for the FRW universe metric, 

3 

ds^ = -dt^ + a{tf {dx') ' , (13) 

i=l 

taking into account the contribution from the quantum anomaly, Ta, the equation of motion ([4]) becomes 

= -^H^ + ^4>^ + V{c^) - 240/'(0)i73 + . (14) 

An equation of this sort (dH) was also obtained for dilatonic gravity [H, [39|. On the other hand, Eq. ^ becomes 

= -7 + 3i?0) - V'{4)) + 24/'(0) [hH^ + H^^ . (15) 

In the above equations it has been assumed that 4> only depends on time. 

Postulating the de Sitter solution, i.e. looking for solutions of the form a — apcxp {Ht), where H and the scalars 
are constant, we have for pA that 

PA^-6b'H. (16) 

After this, Eq. (fT4|) becomes 

= --^H^ + V((j))-6b'H . (17) 

And, since we are looking for solutions with H = const, we are now interested in solving the equations of motion in 
the form 

-V'{^)+2Af{(j)){HH^ + H^) = 0, (18) 
-^H^ + V((t>)-6b'H = 0. (19) 

We may now consider the case when V{(j)) and f{(j)) are given as exponents, with some constant parameters 
Vq , /oand (j)o, 

V = Vo e-^ , /((/)) = /o , (20) 

these are string-inspired values for the potentials. Assuming the scale factor to behave as a = ao exp (Ht), we get the 
following equations of motion 

— l^oe"^ +48 4^ e^i?" = 0, (21) 
9o \(PoJ 

-^H^ + Voe~^ -6b'H = 0. (22) 



Now, combining Eqs. pT|) and ((22|) we find the following equation for H: 



b' 



+ I I ^ + — e ^ = . (23) 



If we switch off the quantum effects {b' — 0), then we find for H^: 

This is the pure classical solution, which was obtained in [23| . This solution may require the coupling /o to be 
negative, otherwise the solution may not exist. Furthermore, since the Hubble rate can be determined by the initial 
condition {H — Hq), the choice of 0o is fully arbitrary. 
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Returning to Eq. (j23p . in looking for a solution which takes into account quantum effects, we easily find a numerical 
one by rewriting Eq. in terms of the potential V, as follows, 

^'^ + + = . (25) 

SfoVoK^ 4Vb/o 

This equation has the form of a reduced cubic equation and thus it is possible to find a solution by using Cartan's 
formula. The discriminant of the cubic equation is 

^=(24)3/3^3,0 (^1 + V j- 

As seen from ((26)) . for the cubic equation (|23p the sign of the discriminant depends on the sign of /o, since we define 
Vo to be positive. 

For the analysis of the solutions of Eq. (|25p in this case, we will better consider numerical solutions and look for 
the behavior of the Hubble rate iJ as a function of the potential V . Let us consider, as a first case, the solution with 
the following values for the parameters: /o = —0.5, Vq = 1, 6' = —0.5, and n = IGttG, where G = 6.63 x 10~^ is the 
gravitational constant. For those, the behavior of H is illustrated in Fig. 1. We immediately see from this figure that 




FIG. 1; Behavior of the Hubble rate H as & function of the potential V, for several values of V. 



we indeed have three possible ways to define if as a function of the potential V. That is, of course, a consequence 
of the three possible solutions of the cubic equation ([25|l . We will be interested only in the case when H has positive 
values for positive values of the potential V. 

As seen in the plot, the Hubble rate H is growing as the potential grows, starting from some minimal positive value. 
From another side, it is important to remark that there is a positive solution for H which is very close to the classical 
solution (dashed line), up to the point where the continuous-line solution tends asymptotically to a very small value. 
This case clearly shows that quantum effects are only a perturbation to the classical solution. It is also clear from 
Fig. [T]the interesting feature that the asymptotic value of H depends on the quantum correction term 6', which on 
its turn depends on the number of fields N. This means that the asymptotic value attains a greater positive value if 
A'' is larger. In the absence of the b' term this discontinuous behavior of the solution for H disappears. This leads to 
the conclusion that the instability of the solution around its asymptotic value is exclusively due to quantum effects. 

It is remarkable that there exists, owing to the quantum effects, another possible solution for H with a positive 
value for /q. This solution exists in a very small region, as a pure effect of quantum contributions since it turns out 
that in the classical case [2^ it is not possible to obtain any positive solutions for H — as seen from Eq. p4p . In Fig. 2 
the situation we have just discussed is shown, corresponding to a positive value for /o. 
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FIG. 2: Behavior of the Hubble rate H for several values of V , with a positive coupling of the scalar field with the Gauss-Bonnet 
invariant; /o > 0. 



III. CONCLUSIONS 

To summarize, a number of quite interesting conclusions can be drawn from the present study of the influence of a 
combination of quantum effects and modified gravity as a possible way for interpreting dark energy in an accelerated 
inflationary de Sitter universe. In particular, it is expected that for large values of the curvature {R — 6H + 12H^) 
the de Sitter epoch thus described corresponds to early-time inflation. 



H 



1x10' 
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1x10"" -2x10 



2x10"" 4x10" 



FIG. 3: Behavior of the Hubble rate H near its value at the dark energy stage. The dashed line corresponds to the positive 
coupling case, /o > 0, and the continuous line to the negative case, /o < 0. 



From another point of view, for small values of curvature [i? ~ (10 eV) ] we expect that the de Sitter era realized 
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in this model will correspond to the dark energy stage. However, as seen in Fig. [3l this occurs for both positive and 
negative values of the potential, depending on the sign of /q. This result clearly shows that, contrary to what happens 
in the classical case, quantum effects lead in fact to the realization of the de Sitter era corresponding to the dark 
energy stage, even for positive values of the /o constant. Such situation may potentially have interesting cosmological 
consequences. 

In a similar fashion, the coupling of GB dark energy with F{R) gravity could be considered, too. This should 
not represent a big problem, in principle, owing to the fact that the curvature is constant on the solutions. That 
possibility deserves further investigation. 
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